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Abstract. The Ha luminosity density of galaxies is used 
to calculate the local diffuse radiation field at the Lyman 
limit as resulting from star formation in the universe. The 
fraction of Lyman continuum radiation leaking out from 
galaxies through their neutral hydrogen is the most un- 
certain parameter of the calculation. A comparison of the 
diffuse radiation predicted from galaxies with that mea- 
sured from all sources of ionization shows that the average 
Lyman continuum escape fraction should be lower than 
1%. This number is lower than the upper limits reported 
so far in a few objects. It is also shown that the luminosity 
density at the Lyman limit resulting from star formation 
is consistent with the luminosity density and the diffuse 
background measured in the far (non-ionizing) ultraviolet. 
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1. Introduction 

It is thought that massive star formation in the universe 
may provide a significant fraction of the background radi- 
ation that maintains the diffuse intergalactic medium and 
the Lyman a forest clouds highly ionized (e.g. Bechtold et 
al. 1987, Songaila et al. 1990, Miralda-Escude & Ostriker 
1990). This contribution would augment that of quasars, 
especially at high redshift where their number density is 
observed to decline. The contribution of quasars itself de- 
pends on how our picture of their evolution is distorted 
by dust obscuration (Fall & Pei 1993). 

Direct observations of the Lyman continuum (Lyc) ra- 
diation escaping from galaxies is however extremely difh- 
cult and, so far, only upper limits have been obtained with 
the Hopkins Ultraviolet Telescope (HUT) in four nearby 
star-forming galaxies (Leitherer et al. 1995). Attempts to 
understand how Ly^ radiation leaks out from sites of star 
formation, ionizes the diffuse interstellar medium around 
and eventually escapes from a galaxy (e.g. Dove & ShuU 



1994, Patel & Wilson 1995a,b, Ferguson et al. 1996) have 
shown that the phenomenon is dominated by patchiness 
in the distribution of the neutral gas and should be highly 
random. Any quantitative assessment of the contribution 
of galaxies to the ionizing background would therefore re- 
quire a large number of observations before a Lyc lumi- 
nosity function is established. This uncertainty on the Lyc 
escape fraction is also a severe limitation for model predic- 
tions even though some of them are reasonably successful 
in linking the Lyc radiation produced by star formation 
to the rate of chemical enrichment in the universe (Cowie 
1988, Songaila et al. 1990, Madau & ShuU 1996). 

In this paper we use the recent determination of the 
Ha luminosity density of nearby galaxies by Gallego et 
al. (1995) to estimate the contribution of galaxies to the 
diffuse radiation background at the Lyman limit and at 
z = 0. As a support to our approach, the Ha luminosity 
density will be compared with other tracers of the local 
star formation activity such as the luminosity density and 
the diffuse background in the far non-ionizing ultraviolet. 
From the comparison between the diffuse radiation at the 
Lyman limit predicted from galaxies and that measured 
from all sources of ionization, we will derive an upper limit 
to the effective Lyc escape fraction in the local universe. 

2. Basic formulation 

Following current formulation (e.g. Bechtold et al. 1987, 
Meiksin & Madau 1993) the mean specific intensity J^^, of 
the diffuse radiation field at frequency i/q, as seen by an 
observer at redshift (outside our Galaxy), writes as 



(1 



:e{iy, z) 



expl 



(l + z)2(H-2go2)i/" 



dz 



(1) 
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where e(z^, z) is the proper volume emissivity (expressed in 
erg cm^'^ s^^ Hz^^) at frequency v — + z) and red- 
shift z, and exp[—Teff{i^o,z)] is the mean transmission of 
a clumpy medium averaged over all lines of sight. In our 
application the volume emissivity is the ultraviolet emis- 
sion resulting from the star formation activity in galaxies. 
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Taking an average over all morphological types and all star 

formation histories, the luminosity density can be consid- 
ered continuous with a unique mean spectral shape. In 
these conditions and following the notations of Bechtold 
et al. (1987); e(i^, z) can be split into e(fo) the current 
local luminosity density of galaxies at frequency vq and 
k{u/u(i) the spectral shape normalized to 1 at the fre- 
quency with a factor V'(^) accounting for any proper 
evolution in the luminosity density (the density variation 
due to expansion is taken into account by a factor (1+2;)^). 
Adopting g'o = 0.5, we get 



i/'(z)fc(— ) , 
J^o (1 + z)^'^ 



exp[-Te//(t'o,2;)]dz (2) 



The effective optical depth Te//(i'o, z) due to Lyman con- 
tinuum absorption of Hi and Heii by discrete absorption 
systems is given by (e.g. Paresce et al. 1980, M0ller & 
Jakobsen 1990, Miralda-Escude & Ostriker 1990) 



(3) 



Teffiuo, z) = dz' / f(N, z') (1 - exp(-r)) dN 
Jo Jo 

where f{N, z') is the redshift and column density distribu- 
tion of absorbers along the line of sight and r the optical 
depth through an individual cloud of column density N. 



3. Stellar contribution to the radiation field at the 
hydrogen Lyman edge 

We now use equation (2) to calculate the radiation back- 
ground at the Lyman limit Jgi2- The luminosity density 
is assimied to have a proper evolution parameterized as 
V-'(z) = (1 + z)'^ . Given the severe uncertainties about 
the spectral energy distribution of Ly^ photons from star 
formation, the normalized spectral shape is modelled as 
k{v/vQ) — {v/vq)'^ = (1 + z)" while we add a factor for 
the absorption of stellar Ly^ photons by the neutral hy- 
drogen of each galaxy. In order to account for the 
freqTiency dependence of the absorption of Ly^ photons 
by neutral hydrogen, we have written the absorption fac- 
tor as ea;}j[-6.3 x 10-'^^Nh{1 + z)-% The equivalent Hi 
column density term {Nh) is introduced for calculation 
purpose (we make no assumption as to the geometry and 
distribution of the gas) and the resulting effective Lyc es- 
cape fraction will be defined as / = ea:;p[— 6.3 x \Q~^^Nh\- 
In these conditions the background radiation J912 (in 
erg cm~^ s~^ Hz~^ sr~^) at the Lyman limit writes as 



Jqi2 = 



(l + z)^'+"-^/2 X 



g-6.3xlO-i«NH(l-Hz)-3 g-reff(912,z) 



(4) 



3.1. Adopted parameters 



By comparison with existing models for the Lyc spec- 
tral energy distribution of star-forming population (e.g. 



Bruzual & Chariot 1993), a value of a in the range — 1 

to —3 is a realistic approximation. An upper bound to 
2; = 2 is also realistic since it cuts any stellar flux con- 
tribution below 304 A . The local luminosity density of 
galaxies e(912) at the Lyman edge can be derived from 
the total Ha luminosity per unit volume of 1.26 x 10"^^ ergs 
s~^ Mpc~^ evaluated for star-forming galaxies in the local 
universe by Gallcgo et al. (1995). Under the current con- 
ditions valid in the ionized gas of galaxies (T=10^K and 
case B) we get a local density of Lyc photons of 9.2x10^'^ 
s"^ Mpc-3 (Osterbrock 1989). This number can be con- 
sidered as a lower limit since Lyc photons in optically 
thin gas produce fewer Ha photons than those in optically 
thick gas. The relation between the Lyc photons density 
and the luminosity density at 912 A depends on the value 
of a. For our average case a = —2, we find a luminos- 
ity density at 912 A of 4.8 x lO^^^ erg s'^ A^^ Mpc^^ qj. 
1.3 X 10^^ erg s~^ Hz~^ Mpc~^. Incidentally, the relation 
logNx/Lx = 13.28 ± 0.16 (photons A erg-^) established 
by Leitherer et al. (1995) would give the same value. This 
relation was established for starbursts with different star 
formation histories and initial mass functions while the 
simplifying assumption of a continuous star formation rate 
is probably valid at the scale of the local universe. 

Significant evolution of galaxies is now well established 
(e.g. Elhs et al. 1996, Lilly et al. 1996, Fall et al. 1996) and 
we adopt 7 = 4 from z = to z = 1 as found by Lilly et al. 
(1996) for the evolution of the luminosity density of the 
universe at 2800 Alnsofar as the light at this latter UV 
wavelength is essentially tracing on-going star formation, 
we think that the same exponent should be valid at our 
shorter wavelengths. Beyond z = 1 the evolution is known 
to slow down but the situation is less certain. We have 
adopted 7 = 0, bearing in mind that this choice is not 
critical since the contribution to the background at z = 
from objects at high redshifts is small as soon as the evo- 
lution is not strong. Last, the calculation is independent 
of the value of Hq since luminosity densities scale as Hq. 

3.2. The intergalactic opacity term 

At the Lyman edge and ignoring absorption due to Hell 
for a line of sight limited to z = 2 (Hei absorption is negli- 
gible), the effective optical depth Te//(912, z) in equation 
(4) writes as 

/•Z POO 

Tef f {912, z) = / dz' / f(N,z') X 
Jo Jo 

(1 - exp [-6.3 X 10-^^(1 + z')-^N]) dN 

Assuming a power-law of exponent —1.5 for the column 
density distribution (Petitjean et al. 1993, Songaila et al. 
1995), and adopting the line densities per unit redshift 
and the evolution parameters from Boksenberg (1995), 
f{N,z') writes as 1.07xlO^N-'^-^{l+z')°-^^ for the Lyman 
a forest clouds, and 5 x 10'^N~'^-^{1 + z')^-^ for the Lyman 
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limit systems {N > 1.6 x 10^'^ cm~^). The calculation of 

the first normalization constant accoimts for the detection 
limit of 0.24 A rest equivalent width and a velocity width 
of 30 km s-i as in Miralda-Escude & Ostriker (1990). The 
parameters for the Lyman a forest clouds have been ob- 
tained for < z < 1.3 but the plot of their evolution up 
to 2: = 3.7 (Boksenberg 1995) shows that our parameteri- 
zation remains appropriate till the adopted limit at z = 2. 
As a numerical example, we find a transmission from z = 
to z = 2 exp{-Teff {912,2)) = 0.384. Playing with the er- 
ror bars given on the line densities per unit redshift and 
the evolution parameters (Boksenberg 1995) we find that 
this transmission does not change by more than 40%. 

3.3. Results 

The background radiation calculated by equation (4) with 
the average transmission models discussed above is dis- 
played in Figure 1 as a function of the Lyc escape fraction. 
Although our evaluation is based on the measured Lyc 
photons density in the local universe and avoids therefore 
most of the uncertainties inherent to pure model calcula- 
tions, it still depends on a few parameters, the index a 
of the average spectral shape in the Lyman continuum, 
the evolution factor and the opacity of the intergalac- 
tic medium. The resulting uncertainties are illustrated in 
Fig. 1. First, the impact of the ill-known index a (values 
—1, —2, —3, —4 are used in Fig. 1), is found to be reduced 
by the relation between a and the luminosity density at 
the Lyman edge for a given Lyc photon density. Second, 
the effect of a larger intergalactic opacity as obtained with 
the upper limits given by Boksenberg (1995) on the den- 
sity of Lyman clouds and Lyman limit systems per unit 
redshift is modest and comparable with a change of one 
unit of the index a. In contrast, the calculation is sensi- 
tive to the amount of evolution as shown with the case 
of a milder evolution 7 = 2 (till z = 2) plotted in Fig. 1. 
Selected as the variable against which the diffuse radia- 
tion has been plotted in Fig. 1, the Lyc escape fraction is, 
as anticipated, the major source of uncertainty. We note, 
however, that the diffuse radiation does not decrease as 
fast as the Lyc escape fraction. The main reason is that 
galaxies at high redshift contribute to the diffuse radiation 
and may still be optically thin at 912/(l + z) A while their 
nearby counterparts are optically thick at 912 A. The is- 
sue of the uncertainty on the luminosity density itself at 
the Lyman edge will be addressed in the two following 
sections. 

4. CompEirison with the ionizing background: con- 
straints on the Lyc escape fraction 

We have compared the predictions of Figure 1 with the 
measurements of the ionizing background (at z f« 0) at 
the Lyman edge. These measurements are only indirect 
and include the contribution of all ionizing sources, chief 
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Fig. 1. Background radiation due to star formation in galaxies 
predicted at the Lyman limit and a = as a function of the 
Lyc escape fraction. Solid lines: spectral index (y'^ , see text) 
a = —1, —2, —3 respectively from top to bottom. Long-dashed 
line: a = —4. Dot-dashed line: a = —2 and intergalactic 
opacity higher than nominal with a density per unit redshift 
(dN/dz)o = 30.9 instead of 24.3 for the Ly forest clouds and 
a density per unit redshift (dN/dz)o = 0.42 instead of 0.25 
for the Lyman limit systems. Short-dashed line: a = —2 and 
(1 + 2)^ evolution from 2; = to 2: = 2 instead of (1 -|- 2)"' 
from 2 = to 2 = 1 and no evolution from 2 = 1 to 2 = 2 in 
all other cases. The horizontal thick line is the upper limit to 
current measurements of the diffuse background radiation at 
the Lyman limit and 2 = 0. 

among them the quasars. They are therefore upper limits, 
possibly very generous, to the contribution of galaxies. 
The measurements have been obtained so far through a 
variety of methods (Bechtold 1993, for a review), such as 
the Ha emission from 21-cm emitting clouds, the mod- 
eling of the sharp edges of Hi clouds and, the proxim- 
ity effect (Kulkarni & Fall 1993), following the obser- 
vations of the Lyman a forest at low redshift with the 
HST. From the recent works of Dove & Shull (1994), 
Vogel ct al. (1995) and Donahue et al. (1995) and re- 
view of previous determinations therein, we conclude that 
the ionizing background at 900 A should lie in the range 
0.06 — 1 X 10^^^ erg cm s^^ Hz^^ sr~-^. At its lower 
bound the ionizing background is comparable to the evalu- 
ations currently made for the contribution of quasars (e.g. 
Madau 1992), indicating a possible negligible role of the 
galaxies. 

Fig. 1 shows that, imder most of the conditions, the 
Lyc escape fraction is smaller than 1%. This is smaller 
than the upper limits in the four star-forming galaxies 
observed by Leitherer et al. (1995), especially after the 
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modifications advocated by Hurwitz et al. (1997) on the 

basis of unaccounted absorption by Hi gas in our Galaxy. 
Our upper limit would be less restrictive than 1% if the 
spectral index a is lower than —4, or the evolution milder 
(7 < 4) than found by Lilly et al. (1996). Insofar as our 
comparison was made with the intergalactic radiation field 
of all origin, our upper limit is probably generous but, even 
if small, does not exclude that galaxies make most of the 
ionizing intergalactic radiation field. 

Our conclusion also depends crucially on the value of 
the total Ha luminosity density of the local universe in 
the sense that the larger the Ha luminosity density (or 
the associated number of Lyc photons), the stronger is 
the upper limit for the Lyc escape fraction. Our limit of 
1% seems firm for two reasons. First, the Ha luminosity 
function built by Gallego et al. (1995) is based on star- 
forming galaxies with Ha+N[H] equivalent widths larger 
than 10 A and therefore should more likely lead to an 
underestimation of the total Ha luminosity density. Sec- 
ond, our conversion of the Ha luminosity density into a 
Lyc photon density under case B assumptions provides, as 
previously said, a lower limit on the latter quantity. 

5. Consistency with the non- ionizing far UV radi- 
ation 

As the Ha luminosity density of the local universe is a 
key input in our calculation, we have tried to evaluate 
how this quantity and the associated number of Lyc pho- 
tons compare with other measured quantities tracing the 
star formation activity in the universe such as the far-UV 
luminosity density as determined from surveys of galaxies 
or the far-UV (non ionizing) background as measured by 
in-orbit experiments. 

The former quantity is related to the luminosity den- 
sity emitted at 900 A by star formation activity through 
two parameters, the Lyman break of a pure star-forming 
population and the extinction at far-UV wavelengths due 
to dust mixed with the young stars in each galaxy. The 
far-UV spectral energy distribution (from longward of the 
Lyman break to « 2000 A) of a pure (without dust) and 
continuous star-forming population can be assumed to 
be flat in energy per frequency unit as shown by the 
models of Bruzual & Chariot (1993) and the observa- 
tions of star-forming galaxies with little extinction by 
Calzetti et al. (1994). In the case of a continuous rate 
as we expect for the average star formation in the lo- 
cal universe, the Lyman break factor should be of the 
order of 4 according to Bruzual & Chariot (1993) and 
possibly between 6 and 20 depending on the initial mass 
function for massive stars according to Leitherer & Heck- 
man (1995). With an average far-UV extinction of the 
order of 1 mag as discussed by Deharveng et al. (1994), 
we get a far-UV (observed) to 900 A (emitted) lumi- 
nosity density ratio between 1.6 and 8. For comparison, 
the 1600 A (observed) to 900 A (emitted) flux ratios in 



the four galaxies observed by Leitherer et al. (1995) are 

foimd to be 0.3, 1.5, 1.6 and 6. With the ratio above and 
our local luminosity density e(912) at the Lyman edge 
of 1.3 X 10^^ erg s~^ Hz~^ Mpc~^ we get a far-UV lu- 
minosity density of 2 - 10 xlO^^ erg s"-^ Hz~^ Mpc"'^. 
This is in agreement with the local luminosity den- 
sity of 3x10^^ W Hz-i Mpc-3 (Ho=50 km s'^ Mpc'^) 
evaluated at 2800 A by Lilly et al. (1996) and of 
6x10^^ W Hz-i Mpc-3 (Ho=50 km s"! Mpc"!) evalu- 
ated at 2000 A by Milliard et al. (1997). 

Although the far-ultraviolet background is rich of sev- 
eral possible components, some of them of galactic ori- 
gin, various arguments have established the accumula- 
tion of galaxy light along the line of sight as the domi- 
nant extragalactic contributor (e.g. Bowyer 1991, Jakob- 
sen 1995, for a review) with an intensity in the range 
50 150 photon cm-2 s"! A"i sr-i(or 0.7 2 xlO^^i 
erg cm~^ s~^ Hz~^ sr~^ at say 2000 A). This background 
radiation can be converted back into a local luminosity 
density e(2000) according to equation (2) which, with the 
same assumptions as above but without the neutral gas 
absorption term and the intergalactic opacity term, writes 
as 

J2000 = --s-e(2000) / (1 + ^)T+«-5/2 dz (5) 

The integral upper bound is now 1.2, the value which 
shifts the Lyman break at 2000 A. Assuming 7 = 4 and 
a in the range —1 to —3 (the case a = is observed for 
galaxies with little extinction and is not appropriate for 
the average spectral shape) the upper limit background 
of 2x10"^"'^ erg cm~^ s~^ Hz^^ sr^^ gives a far-UV lumi- 
nosity density in the range 2.6 - 6 xlO^^ W Hz~^ Mpc~"^, 
again in satisfying agreement with the range of values that 
we have derived above. 



6. Conclusion 

We have used the Ha luminosity density of galaxies in the 
local universe for estimating the contribution of galaxies 
to the ionizing background at z = 0. This approach re- 
duces the number of theoretical assumptions entering this 
type of evaluation and is therefore expected to narrow the 
range of predicted values. Of the remaining parameters 
needed for the evaluation, the fraction of Lyc radiation 
leaking out from galaxies is the most uncertain. Prom a 
comparison between the ionizing photons predicted by on- 
going star formation in galaxies and those measured in the 
intergalactic radiation fleld, it is found that the Lyc escape 
fraction is smaller than 1% in the local universe. This av- 
erage upper limit is more restrictive than the upper limits 
reported in four galaxies by Leitherer et al. (1995) and 
re-analyzed by Hurwitz et al. (1997). 
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